7r-kinks in strongly ac driven sine-Gordon systems 
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We demonstrate that 7r-kinks exist in non-parametrically ac driven sine-Gordon systems if the 
ac drive is sufficiently fast. It is found that, at a critical value of the drive amplitude, there are 
two stable and two unstable equilibria in the sine-Gordon phase. The pairwise symmetry of these 
equilibria implies the existence of a one-parameter family of 7r-kink solutions in the reduced system. 
In the dissipative case of the ac driven sine-Gordon systems, corresponding to Josephson junctions, 
the velocity is selected by the balance between the perturbations. The results are derived from a 
perturbation analysis and verified by direct numerical simulations. 

PACS: 03.40. Kf, 52.35.Mw, 02.30. Jr 



Soliton bearing systems are very important for our un- 
derstanding of collective phenomena in many physical 
systems in the one dimensional approximation |l]]. Of- 
ten, such systems are perturbed in one form or another 
and sometimes these perturbations are temporally peri- 
odic |§-f|]. If the perturbations are small, one can ap- 
proximate the dynamics through the adiabatic perturba- 
tion technique |5|,[| , where the integrability of the unper- 
turbed system is used to assume wave profiles for ana- 
lytical perturbation techniques applied to the "nearly Ln- 
tegrable" system. However, are the perturbations large, 
we cannot simply assume the unperturbed wave profiles 
as representing a good approximation to the dynamics 
since the near integrability is lost. 

This problem was previously addressed for strong peri- 
odic perturbations of sine-Gordon systems J^^), where it 
was shown that the direct ac drive can induce a Shapiro 
type locked phase ||[l0| to which a 27r-kink can local- 
ize. Similarly, it has been demonstrated ||TT[-[T3[| that a 
parametrically driven sine-Gordon system can produce 
stable 7r-kink propagation due to the well known effec- 
tive Kapitza potential of a driven pendulum [flif . 

It is thereby well documented that strongly perturbed 
nonlinear soliton bearing systems can retain some forms 
of near integrability in certain windows of the perturba- 
tion parameter space. Previous analyses for the sine- 
Gordon system have been performed using the analo- 
gies to the single pendulum cases, where Shapiro phase- 
locking exists for the direct drive, leading to stable 27r- 
kinks in a rescaled sine-Gordon chain. 

In this paper we demonstrate that n-kinks can propa- 
gate in strongly perturbed, directly driven sine-Gordon 
chains if the perturbation parameters are chosen near 
the values, leading to zero crossing of the relevant Bessel 
functions determining the size of a Shapiro step. The 
analysis is based on the normal form technique that re- 
lies upon a time scale separation between the rapidly 
oscillating driving force and a relatively slow behavior 



of the residual field. This technique was previously ap- 
plied to the parametrically forced sine-Gordon equation 
(SGE) in [jnJUi- We find for nondissipative ac driven 
SGE a one-parameter family of 7r-kink solutions mov- 
ing with any prescribed velocity. In the case of damped 
and driven SGE, we show that the velocity is selected 
and only one of the 7r-kinks survives. We verified the 
obtained results by performing numerical simulations of 
the ac driven SGE. 

First, we consider the equation of motion for a directly 
forced pendulum cf> + sin0 = Mf(u)t) , where <fi is the 
phase of the pendulum, / is a mean-zero periodic func- 
tion, M is constant, t represents a normalized time, and 
the normalized frequency lo is assumed to be sufficiently 
large compared to the natural frequency of the pendu- 
lum, u>o = 1. We shift the emphasis to the oscillating 
reference frame by the transformation 



= 9 + Muj- 2 F(ujt), 



(1) 



where F has zero mean and F"(t) — /(r). We 
then obtain the parametrically forced equation 9 + 
sin (9 + Aluj~ 2 F(ujt)) — , with the Hamiltonian 



H = y - A(u>t) cos(0) 



B(ut) sin(0), 



(2) 



where p is the momentum canonically conjugate to 9, and 
A(ujt) = cos(Muj- 2 F(Lut)), B(wt) = sm(Mu- 2 F(ujt)) . 
Using 2ir/u> periodicity of A, B, we denote {^4} = A — 
(A) ,{B} = B- (B) , where (...) = (1/2tt) f dr . . . 

Following Refs. [ pdp^ ], we apply the normal form tech- 
nique, to move mean-zero terms to a higher order. Let 
the first canonical transformation be defined implicitly 
as 



p = p 1 + d e W 1 (9,pi,t), 



9 + d Pl Wx(6,p u t). (3) 



The transformed Hamiltonian takes the form Hi = H + 
Wit ■ To remove mean-zero rapidly oscillating terms, we 
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choose W\ — uj~ 1 {A}_i cos(0) — uj~ 1 {B}_i sin(#) , where 
{A}_i is a mean-zero antiderivative of {A}. With this 
choice of W\ the transformed Hamiltonian takes the form 



-w~ l pi ({A}_isin6»i + {£}_! cos 6q) 



(4) 



Unlike the original Hamiltonian (g), the transformed 
Hamiltonian Hi contains terms with small (~ u>~ 2 ) pos- 
itive time-dependent coefficients {A} 2 _ 1 , {B} 2 _ 1 , which 
have non-zero averages for any nontrivial choice of M, f, 
and therefore, cannot be removed from the Hamiltonian. 

These terms have essential effect on the system dynam- 
ics when the lower order potential energy terms vanish, 
i.e. (A) — (B) — (all terms ~ uj^ 1 always have zero 
averages). To remove the explicit time dependence from 
the Hamiltonian up to the terms ~ w -2 , we perform a 
series of canonical transformations similar to (|) (see also 
jjjj). To remove the terms ~ i 
formation 



we apply the trans- 



pi = P2+d ei W 2 (9uP2,t), 9 2 = 6 1 +d P2 W 2 (6 1 ,p 2 ,t) (5) 
with 

W 2 = lj- 2 p 2 {{A}^ 2 wa.(6i) - {5}_ 2 cos(fli)). (6) 

After straightforward calculations we obtain the trans- 
formed Hamiltonian 

H 2 = ^-{A) cos9 2 + (B) sin 6 2 

+ ^ sin 2 9 2 + {{B}\) cos 2 6 2 (7) 

+ ({A}_i{.B}_i) sin26> 2 ) + uj- 2 R + 0(uj- 3 ), 

where R turns out to be mean-zero (R) — 0. Finally, 
applying the third transformation with W3 = oj~ 3 R-\ 
and neglecting terms ~ lu~ 3 we obtain 

pi w -2 
H(P, 6) = — - Ccos(8 - 7) - —D cos(26 - S), (8) 

where 

(A) = Ccos( 7 ), <{A} 2 _ X ) ({B} 2 ^) = 2Dcos(5), 

(B) = CBin( 7 ), -<{A}_i{B}_i) = Dsm(5), 



(9) 



and P and O are new canonical variables. 

For C ^ there is only one stable equilibrium = 7 
and one unstable equilibrium = 7r + 7 , for large 
frequencies; this equilibrium corresponds to the usual 
phase-locked Shapiro state known from Josephson junc- 
tions M. However, as C passes through 0, a bifurcation 
occurs and (for C = 0) the system has two stable equilib- 
ria given by = 5/2, tt + 5/2 and two unstable equilibria 
given by 6 = tt/2 + 8/2, 3tt/2 + 5/2 . 



Now we turn to the directly forced SGE 
<j>tt ~ 4>xx + sin0 = Mf(uit). 



(10) 



After applying the transformation ([j]) we obtain the evo- 
lution equation for a new phase 9 on top of a rapidly os- 
cillating background field (written in the canonical form) 



h = P, Pt 



sm(8 + Mu>- 2 F(ujt)). 



(11) 



Invoking the canonical transformations similar to those 
for the directly forced pendulum (see also [Q) and using 
(H), we obtain for the corresponding Hamiltonian 



H = 



+ DC 



^ + H(p,e) + o(u- 3 ) 



dx, 



(12) 



where the error terms 0(u>~ 3 ) contain the derivatives up 
to the second order. For sufficiently large u>, when we can 
neglect these terms, the above Hamiltonian corresponds 
to the double SGE. After retracing the identical transfor- 
mation (jl]), the obtained approximate solutions become 
7r-kinks on top of the rapidly oscillating background field. 
We now consider the damped and driven SGE 



4>tt - 4>xx + siruf) = Mf(ut) - acpt + f) 



(13) 



which is frequently used to describe long Josephson junc- 
tions where (f> is the phase difference between the 
quantum mechanical wave functions of the two super- 
conductors defining the junction, the normalized time t is 
measured relative to the inverse plasma frequency, space 
x is normalized to the Josephson penetration depth, and 
the nonlinear term represents tunneling of superconduct- 
ing Cooper pairs, normalized to the critical current den- 
sity. The perturbations on the right hand side of the 
equation represent, respectively, a normalized ac driv- 
ing current, a dissipative term arising from tunneling of 
quasiparticles, and a normalized dc bias current. 

To obtain an effective equation of the evolution on the 
slow time scale, we apply the above canonical transfor- 
mations. Since the system is no longer hamiltonian, we 
work with equations of motion rather than Hamiltonians. 
We start with a homogeneous transformation to the os- 
cillating reference frame, <fi — 8 + G(t) , analogous to (ffl), 
designed to remove the free oscillatory term. Substitut- 
ing this transformation to ( |i"3| ) and choosing the function 
G so that it solves the equation G + aG — Mf(ujt) , we 
obtain the equations of motion in the canonical form 

9 t =p, p t = e xx -ap + ri-sin(6 + G(ujt)). (14) 

For the particular case of /(r) = sinr, we find 



G(r) = — 



a M 



M 



w or + lo 



or + uj a 



■ sin r. 



(15) 



Using the notations A = cos G(uit) , B = smG(uit) and 
assuming that (A) = (B) = 0, we apply the series of 
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transformations, as we did for the directly forced pen- 
dulum. This moves all mean-zero terms to higher order, 
leading to 

e t = P + 0(uj- 3 ) (16) 
Pt = ®xx -aP + rj- lj- 2 D sin (26 -S)+0 (w~ 3 ) , 

where D and 8 are given by (|^). Constants a and r\ in 
( pf ) are assumed to be sufficiently small, so that the cor- 
responding terms could be considered as a perturbation. 
Then, in zeroth order in a, 77, the system ( |l6[ ) reduces (af- 
ter neglecting terms ~ ui~ 3 ) to SGE, which has 7r-kink 
solutions. Therefore, slightly perturbed 7r-kinks are ap- 
proximate solutions of the original equation ( |l3| ) on top 
of the rapidly oscillating background field. 

To verify the above predictions, we have performed 
a set of numerical simulations. Rescaling the variables 
in equations of motion ( jiff ) (which in Hamiltonian case 
correspond to the Hamiltonian jl^)), by 



X = 28 - 5, 
x, 



X — s ^ 2 ® 




T= v*± Lt 

UJ 

we obtain the sine-Gordon system 



(17) 



Xt — Pxj Px T ~ Xxx 



In the Hamiltonian case (a = r\ = 0), ( pL8| ) has soli- 
tary wave solutions which, after retracing the transfor- 
mation (p7|), read 6 = V(x, i), P = Vt(x, t), where 



V(x,t) 



2 arctan 



exp 




(19) 



Using (g) we return to the original variables which gives 
us the approximate solution 



= V, p = V t 



{A}. 



- sin(V) 



{B}- 



-cos(V), (20) 



which we use to generate initial conditions for (|ll|). 

We have performed numerical simulations of (11), 
which is equivalent to the original directly forced SGE 
(|l0|). We used a second-order leap-frog method, with ini- 
tial conditions given by ( P0|). The coefficients D and 5 
in (U) are calculated frorn(|): D = 0.2270596 . . . ; i5« 
ir . Figure [l] shows the results of the simulations. The 
driving amplitude M was chosen so as to make both 
(A) and (B) vanish: Muj- 2 w 2.4048 . . .. The veloc- 
ity c in (^) was taken as c = 0.5. The small parame- 
ter used in our perturbation analysis is e = lo^ 1 . One 
can see from the Figure that the kink indeed moves with 
the velocity c 0.5 and that it is stable for t < 300 , 
much longer than e _1 (in the considered case e = 0.1), 
which is a natural estimate for the validity of averaging 
for time-periodic perturbation based on multiple scale 



procedure. For longer times (t > 300), the 7r-kink be- 
comes unstable and eventually disintegrates due to the 
parametric resonances. A similar destruction of the ir- 
kink after a long time was observed in a parametrically 
excited SGE (TIGI 
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FIG. 1. The behavior of 7r-kink solution in the ac driven 
SGE (|Tl]). The solid line corresponds to the initial profile 
given by (|l9|), with c = 0.5 and F{ut) = sin(wt). The dashed 
line corresponds to t — 300. The parameters are ui = 10, time 
step dt = 0.01, mesh size dx = 0.1, system size L = 1500. 



In the dissipative case (a 7^ and r\ ^ 0), we have 



shown that equation (13) reduces to (16), after averaging 
over the fast time scale and neglecting terms of order u>~ 3 . 
Then only one 7r-kink solution of ([if]), with a certain 
value of c, is selected out of the entire family, because of 
the energy balance consideration. 

To find the selected c, we substitute a traveling wave 
ansatz, Z — X — cT, into (p"8|). In the zeroth order in a 
and 77, we obtain \a = 4 arctan [exp (Z/y/1 — c 2 )] . The 
solvability condition for the linearized equation for the 
first order correction Xi> gives us the velocity 



V2lot] f_™ x' dZ 



Ttr/LU 



ay/D J+™{ X ' ) 2 dZ ^/8a 2 D + irhfu 2 



(21) 



We have simulated (M), equivalent to the original 
forced and damped SGE (|l3|), with the initial condi- 
tions obtained from ( p0| ) by adding a small correction 
arcsin(o; 2 ?7/I?) to the first equation of (pfj|), to compen- 
sate an additional energy transferred to the kink from 
the constant source in ( J14[) . The results of the simula- 
tions are given in Figure |2j. The 7r-kink in the Figure 
moves with the velocity c « 0.79 , which approximately 
coincides with the selected velocity c w 0.81 obtained 
from (pl{). The plotted solution remained intact for times 
t < 150. Shortly after that it is destroyed by higher order 
effects such as radiation and nonlinear resonances. 

We have demonstrated that directly strongly ac driven 
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SGE can produce localized 7r-kinks in special regions of 
parameter space. These regions coincide with the regions 
where 27r-kink localization (Shapiro steps) vanishes. The 
formalism for demonstrating the existence of the 7r-kinks 
is based on a time scale separation technique, where the 
ac drive is assumed to be fast compared to any natural 
oscillation in the unperturbed SGE. As a consequence, 
the predicted localization mechanism is limited to the 
high frequency region of parameter space. As the driv- 
ing frequency is lowered, nonlinear mixing between the 
supposed high frequency drive and low frequency wave 
dynamics become more dominant, eventually leading to 
the destruction of localization and all coherent behavior. 
We have numerically tested that long time propagation 
of 7r-kinks is possible for relatively low driving frequen- 
cies, as small as w ~ 5, and moderate driving ampli- 
tudes, M ~ 60. While technical Josephson applications 
of 7r-kink localization seem far removed from the present, 
the above perturbation parameters suggest that the pre- 
dicted effect may exist within experimental parameters 
used for generating Shapiro steps in current- voltage char- 
acteristics of Josephson junctions. 
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FIG. 2. The behavior of 7r-kink solution in the damped 
and driven SGE @. The solid line corresponds to the initial 
profile 6(x, 0) = V(x,0) + arcsin(w 2 r\ / D) , with V given by 
( [l9| ) and G(cut) given by (|l5|). The dashed line corresponds to 
t = 150. The parameters are ui — 6, time step dt — 0.01, mesh 
size dx = 0.1, system size L = 1500, a = 0.03, 77 = —0.003. 
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